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Abstract 

This work puts forward a generalization of the well-known rocking Markovian Brownian ratchets to the realm of 
antipersistent non-Markovian subdiffusion in viscoelastic media. A periodically forced subdiffusion in a parity -broken 
ratchet potential is considered within the non-Markovian generalized Langevin equation (GLE) description with a 
power-law memory kernel 77(f) oc r" {Q < a < 1). It is shown that subdiffusive rectification currents, defined through 
the mean displacement and subvelocity v„, {6x(t)) ~ v^f /r{l + a), emerge asymptotically due to the breaking of the 
detailed balance symmetry by driving. The asymptotic exponent is a, the same as for free subdiffusion, {6x^(t)) oc f. 
However, a transient to this regime with some time-dependent aeffCO gradually decaying in time, a < ass(t) < 1, 
can be very slow depending on the barrier height and the driving field strength. In striking contrast to its normal 
diffusion counterpart, the anomalous rectification current is absent asymptotically in the limit of adiabatic driving 
with frequency Q — > 0, displaying a resonance like dependence on the driving frequency. However, an anomalous 
current inversion occurs for a sufficiently fast driving, like in the normal diffusion case. In the lowest order of the 
driving field, such a rectification current presents a quadratic response effect. Beyond perturbation regime it exhibits 
a broad maximum versus the driving field strength. Moreover, anomalous current exhibits a maximum versus the 
potential amplitude. 
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1. Introduction 

The phenomenon of directed current in unbiased on 
average periodic potentials, such as one in Fig. [T] due 
to a violation of the thermal detailed balance symmetry, 
- the so-called ratchet effect - produced a huge 
literature (see, e.g. Refs. |7, 8] for recent reviews and 
further references). It was studied at length, inspired 
in particular by some biological molecular motors (e.g. 
one-headed kinesins C351) [9-12]. At the same time, 
the fascinating topic of anomalous diffusion lfl3-ll7l 
and, in particular, subdiffusion, where the position vari- 
ance grows sublinearly in time, i.e. {6x^{t)} ~ f with 
< < 1, gained dramatically in the importance and 
interest. This interest is promoted not in the last line 
by the emerging evidence for subdiffusion in such com- 
plex media as actin filament networks, other protein 
solutions, interiors of biological cells 1 18-33], as well 
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as in the conformational dynamics of protein macro- 
molecules 113 1143811 . At the moment, subdiffusion is ex- 
perimentally firmly established in different media with 
different a's in the range of 0.25 - 0.9. It is thus natu- 
ral that these two paradigms-breaking(-and-forming) re- 
search lines are currently crossing. 

Whether thermal ratchets based on anomalous subd- 
iffusion are possible or not presents a highly nontriv- 
ial challenge, even on the level of basic and some- 
what oversimplified models. In particular, a periodically 
forced continuous time random walk (CTRW) subdif- 
fusion and the associated fractional Fokker-Planck dy- 
namics ll39l - l42ll turned out to be not sensitive to the ac- 
tion of time-periodic fields in the long time limit with 
an exception of enhancement of unbiased subdiffusion 
1 41 - 43I1 • This fact seems to rule out the very possibiUty 
of rocking (i.e. fluctuating tilt |[lll2l) subdiffusive ratch- 
ets based on the CTRW mechanism Hi- 15, 44, .45J, and 
possibly also on the quenched disorder llT3il46ll . leaving, 
however, a door for the flashing (i.e. pulsating potential 
lA^) subdiffusive ratchets based on these mechanisms 
open ll47ll . 
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In biological applications, the subdiffusion scenario 
based on the alternative to CTRW, fractional Brown- 
ian motion (FBM) [48. 49.1 and associated generalized 
Langevin equation (GLE) descri ptio ns ll36 . can. 



however, be more relevant 11291 15411 . in particular, be 



cause of such a subdiffusion is ergodic ll54i , .55.1 . It does 
not require a quasi-infinite mean residence time in any 
tiny spatial domain, or trap [54], and in contrast with the 
CTRW subdiffusion [56] does not suffer from such non- 
ergodic features as spontaneous immobilization of some 
particles for the time of observation [45]. Given a small 
number of functionally specific (sub)diffusing macro- 
molecules, a sudden "non-ergodic" standstill could have 
fatal consequences for the cell functioning. 

However, what is the physical origin of GLE sub- 
diffusion? In complex environments, such as interior 
of biological cells densely stuffed with different macro- 
molecules, including actin filaments, either mobile, or 
building up the cytoskeleton networks, static, as in eu- 
caryots, or dynamically flickering, as in bacteria, it can 
be due to long range negative correlations in the dif- 
fusing particle displacements (and velocity alternations) 
caused by crowding. A particle moving with an in- 
stant velocity v(f) in certain direction invokes besides 
a viscous drag also a quasi-elastic response of environ- 
ment opposing the motion. The total dissipative force 
can be written in linear approximation as F^-eiit) - 
- j^^ rjit - t')v{t')dt' with a memory friction kernel r]{t). 
For an oscillatory motion, v(f) ~ via)) exp(-/wf) + c.c, 
with frequency a>, the friction force becomes frequency- 
dependent, with amplitude Fy^eii^^) - 
where fjio)) = J^^ rjit) exp(ia>t)dt is the frequency- 
dependent friction Q. The subdiffusive behavior can 
emerge for some special forms of //(f) (see below) when 
this friction diverges at zero-frequency, i.e. ^(0) - 

poo 

Jq ri(t)dt — > oo, or also on some sufficiently long tran- 
sient time scale, when 77(0) is large but finite, which is a 
more realistic assumption. Furthermore, the fluctuation- 
dissipation theorem (FDT) dictates that at thermal equi- 
librium the corresponding energy loss due to dissipation 
must always be compensated on average by an energy 
gain due to unbiased (on average) random force of envi- 
ronment ^(t). For a particle of mass m initially localized 
(y(t ) = for f < to = 0) this leads to GLE description 
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mx + I Tjit — t')x(t')dt' + 



dV(x, t) 
dx 



'it is proportional to a frequency-dependent viscosity, ({(li), 
fi(Li>) cx: l{co), which was introduced in the theory of viscoelasticity 
by A. Gemant fSTll along with the formalism of fractional derivatives. 



where the memory kernel and the autocorrelation func- 
tion of noise are related by the fluctuation-dissipation 
relation (FDR) \M 



mm)^kBTrii\t-t'\) 
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at temperature T . The thermal random force ^(f) has 
necessarily Gaussian statistics within the linear friction 
approximation considered here 11641]. i. e. within linear 
FDT, but not necessarily otherwise I 65ll . 

In the following, V{x, t) - U{x) - xf(t), includes 
some periodic spatial potential U{x) - U{x+L) with pe- 
riod L. It is modulated in time by an external force /(?)• 
In the numerical simulations below it will be considered 
harmonic, fit) - A cos(Qf), with amplitude A and angu- 
lar frequency Q. Apart from above phenomenological 
justification, the GLE description can be derived from 
the Hamiltonian dynamics ifssl IHill . i.e. from the first 
principles. All this makes the GLE approach to anoma- 
lous diffusion and relaxation processes ever more attrac- 
tive. For spatially and temporally unbiased (on average) 
force fields f{x,t) - -dV(x,t)/dx, a rectification cur- 
rent (i.e. ratchet effect) can emerge if only the symmetry 
of detailed balance is destroyed by an external driving 
^ which provides also an energy supply to drive the 
rectified dissipative motion in certain direction. 

A popular model of viscoelasticity with 



77(f) = /7„r7r(i-a), 



(3) 



where Q < a < I (the factor r(l - a) is to relate 
our description with other^, was introduced by Ge- 
mant 0[57[. It allows to rationalize the Cole-Cole di- 
electric response of harmonically bound, U(x) - kx^/2, 
overdamped particles in the inertialess limit, m — > 
11681 1691 . Moreover, this model corresponds to the so- 
called sub-Ohmic spectral bath density 7(w) oc t/^w" in 
the language of the Hamiltonian system-bath descrip- 
tion leading to GLE upon integration over the initially 
canonicaUy distributed bath variables at temperature T. 
The corresponding relation between 77(f) and J(a)) is 
77(f) = (2/n) duJicD) cos(wf)/w 



^Then the frictional term of GLE can be abbreviated as rja ^-^^ 
making use of the definition of fractional Caputo derivative of the frac- 
tional order a |67]. For this particular kernel, the corresponding GLE 
is named sometimes fractional. 

^More precisely, he introduced a more general model for the com- 
plex viscosity which in a particular case yields = f„T'""/[l + 
(-('(xiTc)'""]. The corresponding r](t) oc f(r) matches Eq. ^3) for 
t ^ Tc, while for t » its asymptotics is ri{t) oc f'^^. In Ref i63l . 
Cole and Cole have remarked that the choice l(iS) = fol- 
lowing to Gemant (stIi (our notations are different) leads to their now 
famous form of the dielectrical response function. This corresponds 
to the limit Tc — ► 00. 
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Within this model the variance of free subdifFusion, 
or subdiffusion biased by a constant force evolves in 
time asymptotically (for m — > 0, exactly if to assume 
initial velocities thermally distributed) as 1 5^, 51 ] 



<&2(f))~2D[°'f7r(l+a), 



(4) 



where is the subdifFusion coefficient, which is re- 
lated to the anomalous friction coefficient rj^ and tem- 
perature T by the generalized Einstein-Stokes relation. 

The strict power law kernel represents, however, 
rather a theoretical abstraction. All the realistic power 
law kernels have cutoffs. A particular functional form 
with (upper) incomplete gamma-function. 



7/(0 = aijaT-'Ti-a, f/T,)/r(l - a). 



(5) 



corresponds to the Davidson-Cole dielectrical suscepti- 
bility of the harmonically bound particles ll69ll which is 
also typical for complex fluids, gels and glasses. For 
t « Tc, this kernel coincides with one in Eq. (O. How- 
ever, for f > Tc it has an exponential cutoff which makes 
the zero-frequency friction ^(0) = arjaTlr" finite. In 
such a more realistic case, subdifFusion occurs on the 
time scale Tb < t <^ Tc and turns over into normal dif- 
fusion for t > Tc- Tc can be, however, large enough (e.g. 
seconds to minutes, as in the interior of biological cells) 
for the subdifFusion-limited reactions to become impor- 
tant, in particular because of a finite size of biological 
cells since macromolecules can subdifFuse over the cell 
volume within some time less than Tc- Initially difFusion 
is ballistic on the time scale < f < t;, due to inertial 
efFects, where ballistic time depends on the details of 
memory kernel and mass m. It physically corresponds to 
the relaxation time scale of the momentum (on this time 
scale the ballistic superdiffusion is persistent). Interest- 
ingly enough, a short-time superdiffusion was already 
experimentally observed in viscoelastic fluids |70]. 

The noise ^{t) corresponding to the model in Eq. (O 
is mathematically the fractional Gaussian noise (fGn) 
by Mandelbrot and van Ness (a^. It provides an impor- 
tant instance of the so-called 1 // noises which encom- 
pass noises with a low-frequency power law feature in 
their power spectrum, S{(jS) oc < y < 2 19311 . 

here y - I - a. Moreover, in the absence of po- 
tential the considered model reproduces in the Umit 
m the fractional Brownian motion (FBM) with 
anti-persistence |48, 49J. The motion in periodic poten- 
tials is also antipersistent and subdiffusive js^ll . More- 
over, it has been recently noticed that for a < Qc ~ 0.41 
the inertial effects are generally important ||7l|. With 
their inclusion, this model explains also a generalized 



Rocard dielectric response ll72ll of harmonically bound 
particles (r^ oo) which can account also for the fast 
y6-relaxation due to the cage effect. A corresponding 
generalization of the Davidson-Cole succeptibility to in- 
clude inertial effects can also be readily given (as well 
as for any other memory kernel and model of anomalous 
dielectric relaxation and "glassy" behavior). 

For particular viscoelastic media, the corresponding 
memory kernels 77(f) can be derived from the rheolog- 
ical measurements ifisl. T9I 21, 28]. This is a phe- 
nomenological viewpoint. Moreover, one can derive 
many memory kernels from a theoretical approach, e.g. 
from the mode coupling theory |61] or from the the- 
ory of viscoelasticitity of semiflexible polymer entan- 
gled networks L73.1 . The later one yields e.g. a power 
law kernel with a - 0.75 which agrees with some exper- 
iments (l^Jl^, and possesses also a long-time memory 
cutoff naturally related to the polymer network proper- 
ties. 

Moreover, such a GLE subdifFusion with a = 0.5 
has been used to mimic [741 a potential-free single- 
file subdiffusion of hard core Brownian particles due 
to a geometric crowding effect (i.e., a particle cannot 
move farther than some typical distance until another, 
road-blocking particle moves) ||75-77]. The latter one 
has been experimentally realized and observed with col- 
loidal particles 178., ,79,1 . Even if the presence of a bi- 
asing force acting on all particles clearly destroys this 
correspondence (since all the particles drift normally in 
the same direction), one can imagine a situation where 
an external field acts only on some, e.g., electrically 
charged guest particles. Then one can expect that this 
correspondence might hold further and the behavior of 
the charged tracer particles (e.g. a globular proteins in 
a suspension of neutral colloidal particles confined in a 
quasi-one-dimensional geometry) in some electrical po- 
tential V(x, t) externally imposed can yet be mimicked 
by GLE in agreement with |74, S^j- Similar colloidal 
systems |7^ 79.1 can be considered as plausible can- 
didates for an experimental realization of subdiffusive 
ratchets with a - 0.5 which we consider in the fol- 
lowing for the purpose of illustration. The qualitatively 
same efFects were also found for a - 0.75 and are ex- 
pected for other values of a as well. 

The emergence of subdiffusive ratchet effect is highly 
nontrivial. It has been recently shown numerically for 
sinusoidal potential fsT] that such a viscoelastic sub- 
difFusion is asymptotically not sensitive to the ampli- 
tude of the periodic potential, in agreement with 1181 1 
0. This fact might seem to rule out the very possibil- 



The model in Ref. [81] is quantum-mechanical and it includes 
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Figure 1: Considered ratchet potential, Eq. ^6) witliout bias and under 
the critical forward tilt (potential barriers disappear for the motion to 
the right), as well as under the opposite tilt of the same magnitude. 
Take a note that the tilt in the opposite direction does not result in a 
barrierless total potential. This intuitively selects (think about an over- 
damped particle sliding in the tilted potential) the right direction for 
the net current in the case of a slow time-periodic forcing, symmetric 
and unbiased on average. For a smaller rocking force amplitude, the 
potential banders in the right sliding direction are typically smaller. 



ity of such a rocking ratchet as well. On the contrary, 
we show in this work that such subdifFusive Brown- 
ian ratchets are possible and moreover they display a 
resonance-like dependence on the driving frequency Q.. 
The rectification current vanishes in the limit of adia- 
batic rocking, in a clear contrast with the normal diffu- 
sion rocking ratchets, where it is maximal jTlHtl- This 
behavior somewhat resembles pulsing potential ratch- 
ets of normal diffusion. The physical mechanisms are, 
however, different. In particular, our ratchet also ex- 
hibits the phenomenon of (sub)current inversion at suffi- 
ciently fast driving, similar to its normal diffusion coun- 
terpart 12J, whereas the subdiffusion coefficient remains 
rather robust and weakly sensitive to the details of po- 
tential and driving, being close to that of free subdif- 
fusion. With a further increase of driving frequency 
the rectification current gradually vanishes. We shall 
demonstrate below numerically and explain how this 
puzzling non-adiabatically rocking, anomalous subdif- 
fusive ratchet mechanism operates. 



tunneling effects. In the quantum case, this result is also supported by 
the so-called duality transformation mapping between the dissipative 
tight-binding and full potential problems with different bath spectral 
densities 1621 . However, our results show that this puzzling effect is 
of purely classical origin. It is invoked by the anti-persistent nature of 
GLE subdiffusion. 



2. Theory 

For arbitrary potentials, including a typical ratchet 
potential ||2l 

U{x) = -Uo{^m{2nxlL) + (1/4) sin(47rjt:/L)] (6) 

with period L and height f/o, which we consider hence- 
forth, the exact non-Markovian Fokker-Planck equation 
(NMFPE) which corresponds to the considered GLE 
model is not known. The cases of linear and parabolic 
U (x) present the only known exceptions ll82l487ll which 
cannot much help in the context of ratchet problem. 
This, in particular, makes a rigorous analytical study of 
subdiffusion in nontrivial potentials a highly notrivial 
problem without exact analytical solutions. 

We approach it numerically by embedding the 
non-Markovian subdifFusive dynamics as a Markovian 
stochastic dynamics into the phase-space of higher di- 
mensionality following to Ref. jsill . Namely, the con- 
sidered power law kernel is approximated by a sum of 
A'-exponentials, 



N-\ 



V(f) = ^ rji exp(-v/f), 



(7) 



with V/ = V{)lb' and 7// = {rialT{\-a))Ca{b)v^lb"' scaled 
hierarchically using a dilation parameter b > 1 . In the 
theory of anomalous relaxation similar expansions are 
well-known lll3l 18811 . In the present context, the ap- 



proach corresponds to expansion of fractional Gaussian 
noise into a sum of uncorrected Ornstein-Uhlenbeck 
(OU) noises, ^(f) — YIiJq ^i(t), with autocorrelation 
functions, {(i(t)^j(t')) = kBTriiSije\p(-Vi\t - t'\). Then, 
the GLE ([T]| with the memory kernel in O can be ob- 
tained by eliminating the auxiliary variables m,(0 from 
the following Markovian stochastic dynamics in the 
D - N + 2 dimensional phase space (x, v, mq, un~\)'- 



mv — — - 



d V{x, t) 
dx 



N-\ 



i=0 



iii = -riiV - ViUi + yjlvirjikBT ^i(t) , 



(8) 



where ^,(f) are independent unbiased white Gaussian 
noise sources, {^i(t)^j(t')} = ^^^(f - t'). To enforce the 
FDR ^ for all times, the initial m,(0) are independently 
distributed with the standard deviations cr,- - yjkBTrji 
and zero mean 1 51 , 54 , g^l • The idea of such a Marko- 



vian embedding of non-Markovian GLE dynamics is 
pretty old [89] and the embedding is not unique 151, 90ll. 
However, our scheme is the simplest one which serves 
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the purpose and it leads to an insightful and simple in- 
terpretation of the anomalous rate processes in terms of 
slowly fluctuating rates, with the rms amplitude of rate 
fluctuations which is gradually dying out upon increas- 
ing the potential barriers. For sufficiently high poten- 
tial barriers (e.g. exceeding 12 ksT for a - 0.5 and 
about 9 ksT for a = 0.75 05411 ) the rate description for 
the escape processes is restored. Then the escape rate 
becomes excellently described by the celebrated non- 
Markovian rate flieory ll60[l86ll9lll92ll. 

Moreover, independently of the FBM connection, the 
model in Eq. ([8]) can be considered as a physically plau- 
sible viscoelastic model in its own rights as any power 
law dependence observed in the nature can be approx- 
imated by a finite sum of exponentials and such an ex- 
pansion is a standard methodology e.g. in spectroscopy 
(approximation of spectra by a sum of Lorentzians) and 
in modeling of ion channel kinetics in biology. The in- 
verse of vo corresponds to the fastest time scale of the 
physical noise, or the high-frequency cutoff introduced 
into the (thus approximated) fGn, and vai_i = vq 
corresponds to the low-frequency cutoff. Such or sim- 
ilar cutoffs are always present in any realistic physical 
setup [93]. Even if on the time scale t > b^^^ /vq = Tc 
the subdiffusion will turn over into the normal diffu- 
sion, this time can be extremally large and practically 
not reachable, as in our simulations (or sufficiently large 
for subdiffusion to become physically relevant, as in real 
experiments). As a rule of thumb, a power law extend- 
ing over time-decades can be reasonably well approx- 
imated by a sum of about exponentials choosing the 
dilation parameter b = 10. The number can thus be 
surprisingly small. For example, an experimental power 
law extending over about 5 time decades was nicely fit- 
ted to a sum of 6 exponentials in Ref. f^T] in a more 
flexible way, i.e. not assuming a precise scaling |fl So, 
choosing vo = 10^ (arbitrary units) and = 16, one 
can fit with Ci/2(10) = 1.3 the power law kernel for 
a - 0.5 from t - 10"^ tifl t - 10", i.e. over 14th time 
decades lf54ll . The "error" (if to think of the model of 
fractional Brownian motion with inertia as something 
really fundamental which is not so) introduced by such 
an approximation is within 1% only (5% . Moreover, 
in any realistic experimental setup with some observed 
viscoelastic memory kernel exhibiting normally differ- 
ent complex patterns, see e.g. in Refs. lfl8ll28tl . one can 



Of course, fitting a power law by a sum of exponentials does not 
answer the question about which physics is behind the power law scal- 
ing. However, it oft'ers the way for an economical Markovian embed- 
ding of the observed non-Markovian dynamics. Moreover, the sum 
of exponentials can reflect some real hierarchical, tier like structured 
relaxation with a finite number of tier levels fl^ . 
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Figure 2: (Color online) Stochastic simulations with embedding di- 
mension Z) = 18 for difl'erent potential amplitudes (in units of ksT). 
The initial velocities are thermally distributed and initial diffusion is 
always ballistic. It changes into subdiffusion which is asymptotically 
not sensitive to the barrier height. The transient to this asymptotic 
regime can be however very slow, depending on the barrier height. 
For Uo = 3, it is still not completed, aeffW ~ 0.73 for < t < 10'*. 
It will but decrease to a = 0.5, because of the subdiffusion in peri- 
odic potential cannot be faster than free subdift'usion, even if we can- 
not arrive at this regime numerically within two weeks of simulation. 
Stochastic Heun algorithm is used with time step Af = 10"'', a = 0.5, 
r = 10, and n = 10* trajectories for the ensemble averaging. 



do another expansion, not assuming any precise scal- 
ing, but rather simply fitting the experiment. Our flexi- 
ble methodology will work anyway and the embedding 
dimension D - N + 2 can typically be even smaller 
For example, for all observed cases of biological sub- 
diffusion (typically several time decades) the choice of 
A^ = 6 is sufficient (will be published elsewhere). We 
are interested here but in a faithful Markovian embed- 
ding of the FBM type of subdiffusion (which formally 
requires infinite dimensions) and integrate GLE using 
embedding dimension D = 18 (A^ = 16) which provides 
excellent approximation. 

3. Results 

It is convenient to scale the coordinate x in units of L, 
and the time in units of t = [L^ria/iksT)] so that the 
particle freely subdiffuses over a distance of the order L 
in time t. More precisely, in this scaling D^* - 1 for the 
strict power law memory kernel at any temperature. The 
potential amplitude parameter Uq is scaled in the units 
of ksT and the driving force in the units of ksT/L. The 
frequency is scaled in the units of t"'. Furthermore, the 
influence of the inertia effects is captured by the param- 
eter r = r/Ti, where Tb - L/vj is the ballistic time for 
vanishing friction and vt - y/kBT/m is the rms thermal 
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velocity. High damping corresponds to r » 10. We 
used the stochastic Heun algorithm |95] to integrate the 
system of stochastic differential equations It is of 
the second order of weak convergency in the integration 
time step Af for our particular case of additive noise. 
The Mersenne Twister pseudorandom number generator 
was used to produce the uniformly distributed random 
numbers which were transformed into Gaussian random 
numbers in accordance with the Box-Muller algorithm. 
With the time step At - 10"^ and n - 10* trajectories 
used for the ensemble averaging, the typical accuracy 
of our simulations is within margin of several percents 
ll5411 . as tested by comparison with the exact analytical 
results available for the potential-free subdiffusion and 
for parabolic potentials. This is a very good quality for 
stochastic numerics. It must be also noted that the use 
of double precision floating-point arithmetics cannot be 
avoided to arrive at convergent results. 

3.1. Static ratchet potential 

The influence of static ratchet potential U{x) on the 
subdifussion is illustrated in Fig. |2] After a short 
ballistic stage (within a potential well), the diff'usion 
can look (depending on the potential barrier height) 
initially closer to normal. This is because of a fi- 
nite mean residence time in a potential well exists 
and the escape kinetics, being asymptotically stretched- 
exponential, tends gradually to the normal exponential 
kinetics with an increase of the potential barrier height 
II54II . However, its slows down and asymptotically, in- 
dependently of the barrier height (which is about 2.2Uq 
for static potential), reaches the boardline of free sub- 
diffusion which clearly cannot be crossed. This typical 
behavior can be characterized by a time-dependent ex- 
ponent Q'eff(0 defined from the slope of {6x^{t)) curve 
in the double-logarithmic coordinates: Ceff — > a with 
t ^ 00, independently of Uq. Such a behavior reflects 
the physical nature of viscoelastic subdiffusion which 
is due to the long-ran ge c orrelations in the particle's 
coordinate increments f54j, being in a sharp contrast 
with the semi-Markov CTRW subdiffusion, where such 
correlations are absent in principle. These are not the 
rare events of the escape from potential wells which de- 
termine asymptotically the temporal pace of diffusion. 



'Another useful scaling is: time in r,. = (m/;;,,)'"^ and energy 
in £0 = L^thi^^ The temperature is scaled then as f = 

kgT IEq. The connection between these two scalings is provided by: 
Tv = rr"-^/^-^^"^ and f = r"^""^""'. One can easily rescale results at 
the fixed temperature (like in this paper) from one scaling to another 
one. Other scalings are also possible. 
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Figure 3: (Color online) Mean particle position as function of time for 
difi"erent driving frequencies, f/o = 3, A = 4. Other parameters are 
the same as in Fig. |2] The subvelocity v„ is obtained by fitting the 
numerical dependences by v„?''^/r(l + a) within the last window of 
length AT" = lO' in simulations (as it is shown for Q. = 5). Notice that 
this provides only an estimate since still a^g(t) a: 0.58. The smaller 
is Uq, the more reliable is this estimate, since ffefj is better relaxed to 
the asymptotic value a = 0.5. 



but the anti-persistent nature of viscoelastic subdiffu- 
sion which finally wins and limits diffusion in the pe- 
riodic potentials by the free subdiffusion limit which 
is finally attended regardless the potential height. The 
transient to this astounding asymptotical behavior is, 
however, extremally slow and for this reason may be 
not achieved in practice. Even if the presence of peri- 
odic potential does not influence subdiffusion asymptot- 
ically (and therefore the adiabatically rocking subdiffu- 
sive ratchets are simply impossible), it does profoundly 
influence the whole time course of diffusion. This fact, 
which seems to go beyond any analytical treatment, is 
at the heart of anomalous ratchet effect. This physi- 
cal picture implies that the subdiffusive rectified current 
has a resonance-like dependence against the driving fre- 
quency Q, because of in the hmit Q — > 00 the ratchet 
effect must asymptotically vanish as well. 

3.2. Rocking ratchets 

To verify these qualitative considerations, we have 
tilted periodically the ratchet potential in Eq. ^ forth 
and back by fit) - A cos(Or)- The amplitude A is scaled 
in the units of ksT/L. For a given U[) (in units of hgT) 
there are two critical values of A for maximal force /(f): 
(i) Ai = {3n/2)Uo ~ 4.71 f/o (when the potential barrier 
vanishes for the forward maximal tilt) and (ii) A2 -2A\ 
(same for the backward tilt). Moreover, for A between 
A3 - jtU() and Ai, the potential U(x) has two small 
barriers within each spatial period of V{x) for the for- 
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ward tilt. Generally, for a subcritical driving the po- 
tential barrier towards the direction of tilt is smaller for 
the forward tilt than for the backward tilt of equal am- 
pUtude. Therefore, intuitively the rectification current 
should flow in the positive direction for Uq > 0, cf. Fig. 
[1] and this is indeed the fact, cf. Fig. [3] for a sufficiently 
slow driving. However, for the normal diffusion ratchets 
the current inversion is possible for a sufficiently high 
driving frequency Q |2]. A similar counter-intuitive in- 
version was found also in present case, see in Figs. I3I4I 
Subdiffusive currents we define through the subveloc- 

ityEI 



r(l + a) lim 



{x(t)) 



(9) 



Likewise, the subdiffusion coefficient in the potential is 
defined as 



Da = + a) hm — - — , 

-,(0) : 



(10) 



so that it should coincide with D„ in the absence of po- 
tential. Furthermore, the quality of the anomalous rec- 
tification can be characterized by a generalized Peclet 
number which we define as Pe,y - VaL/Da by analogy 
with the normal diffusion case 19611 . 

It must be noted that practically we are dealing with 
some time-dependent ffeffCO in our numerics which re- 
laxes very slow (depending on Uq) to the asymptotic 
value a. This asymptotic value is not easy to obtain for 
high Uq. To arrive at the end point f = 10^ in our sim- 
ulations for one set of parameters and to have numer- 
ically reliable convergent results one has to propagate 
dynamics for about two weeks (At - 10"'*, n - 10'^) on a 
single node of our modern Linux cluster using a highly 
optimized for the "number crunching" FORTRAN In- 
tel compiler. We derive the corresponding values of Va 
and Da by fitting the dependencies (x(f)) and {Sx^(t)} 
with the afi-^ dependence (extracting the corresponding 
a) within the last time window of length AT - 10^ in 
the simulations. For Uq < 2 this gives reliable results, 
and for Uq = 3 the numerically derived values are less 
reliable (see but a fitting in Fig. |3]for Q = 5 to get an 
idea on the quality of approximation). This is because of 
the corresponding Q'efF(0 is still in fact about 0.58 within 
this time window, if to consider o-eff as an independent 
fitting parameter, i.e. it still did not relax completely 
to a. This value is but essentially lower and closer to 



'The limiting procedure should be understood here in a physical 
sense, i.e. f is large but still much smaller than the memory cutoff 
Tc- Otherwise, the proper mathematical limit will yield infinity for a 
memory kernel with cutoff. 
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Figure 4: (Color onUne) Anomalous current (subvelocity iv) as func- 
tion of the driving frequency for different Uo and A. Rectification 
vanishes for both adiabatically slow and for a very fast driving. A 
cun'ent inversion is observed for sufficiently large CI in all cases. 



a than the corresponding o-eff ~ 0.73 of the non-driven 
subdiffusion in the same potential, cf. Fig. |2] Generally, 
periodic driving accelerates the ultraslow relaxation of 
Q;eff(0 to a. However, for J/q = 3 it is, strictly speaking, 
still not quite achieved. We can compare only asymp- 
totic transport coefficients in different potentials (a com- 
parison makes just a little sense for some snapshot val- 
ues aetf(O)- For larger Uq (e.g. Uo = 4), it becomes 
practically impossible to estimate these quantities. For 
this reason, such cases were not feasible for a quanti- 
tative study, even though the anomalous ratchet effect 
exists of course also. This striking feature of a time- 
dependent Q'efF(0 must be kept in mind when to study 
subdiffusive ratchets experimentally. 

The absence of asymptotic ratchet effect in adia- 
batic driving limit (a transient one yet exists!) is in a 
sharp contrast with the normal diffusion rocking ratch- 
ets where the current is maximized for Q — > 
An inspection of single particle trajectories (not shown) 
indicates that the optimal frequency in our case corre- 
sponds to a synchronization of the potential tilts and the 
particle motion betwe en p otential wells being a kind of 
stochastic resonance n97\\ . In this respect, our anoma- 
lous rocking ratchet resembles more normal flashing 
ratchet than normal rocking ratchet. Indeed, if dur- 
ing a driving (half-)period the particles can move over 
many lattice periods the potential profile asymptotically 
does not matter and the current is suppressed. For the 
frequencies much larger than this inverse characteristic 
time the current is again suppressed (here per an anal- 
ogy with normal diffusion ratchets: too frequent force 
alternations hamper the motion). This explains the oc- 
currence of the resonance like feature. The maximal 
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free subdiffusion limit 
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Figure 5: (Color online) Subdiffusion coefficient as function of driv- 
ing frequency for different Uo and A. 



Figure 6: (Color online) Anomalous current as function of driving 
amplitude. Quadratic response fit is also displayed for comparison. 
J/o = 2, a= 1. 



amplitude of the rectification current first increases with 
increasing the barriers height, see in Fig. |4] where the 
current is maximal for the fixed A at the highest poten- 
tial barrier being optimized for O. It will however be- 
come obviously suppressed with a further increase of 
U(). There are optimal potential amplitudes depending 
on driving strength and frequency, just like in the case 
of normal ratchets. The unexpected analogy to the flash- 
ing normal ratchets is, however, not complete as the 
(sub-)current inversion for sufficiently high frequency 
Q. shows. This feature is similar to one in rocking nor- 
mal ratchets and it unfortunately cannot be explained in 
simple intuitive terms. 

The anomalous diffusion coefficient does not display 
such a profound dependence on frequency as Vq,. It be- 
comes somewhat increased (less than few percent for 
A = 4 and no more than 20% for A = 6) as compare 
with the free subdiffusion coefficient (whose numerical 
value Da ~ 1.012 agrees well with the theoretical value 
of D^* = 1 in the used scaling). However, this effect 
is relatively small, compare Fig. |4] and Fig. |5] This 
remarkable feature is of the same origin as insensitivity 
of the GLE subdiffusion to the height of periodic poten- 
tial. It is in a sharp contrast with the normal diffusion 
case. For this reason, the dependence of the generalized 
Peclet number on frequency just resembles the behavior 
of the absolute value of subcurrent in Fig. |4] 

The dependence of rectification current on the driv- 
ing field amplitude A is shown in Fig. |5]for Uo = 2 
and Q = 1. It displays the same quadratic nonlinear 
dependence on driving amplitude for a small driving 
strength as in normal ratchets Indeed, the rec- 



tification current cannot be present within the linear re- 
sponse (due to the Onsager symmetry of kinetic coeffi- 
cients). It first emerges as the lowest, second order non- 
linear effect which is not forbidden by symmetry con- 
siderations. Moreover, beyond perturbation theory the 
rectification current displays a very broad maximum for 
Ai < A < A2, in accord with intuition. 

4. Discussion and Summary 

In this work we put forward a model of subdiffusive 
Brownian ratchets within the GLE description. Notably, 
a strict power law kernel is not required to reproduce 
subdiffusion on a very long time scale (the theoretical 
asymptotics of normal diffusion for t > 10" cannot 
be even reached numerically on the available comput- 
ers what makes it practically irrelevant in our simula- 
tions), and the Markovian embedding dimension can 
be surprisingly small, allowing to accurately approxi- 
mate nonlinear subdiffusive dynamics driven by Gaus- 
sian 1 noise. Moreover, power law memory ker- 
nels with cutoffs are in fact more physical than strict 
power laws and their expansion over exponentials can 
reflect the corresponding relaxation spectrum of vis- 
coelastic response of the surrounding medium. 

The way we did this Markovian embedding, i.e. ap- 
proximating the fractional Gaussian noise by a sum of 
Ornstein-Uhlenbeck processes is also insightful. In par- 
ticular, it allows to justify the view of the escape dynam- 
ics out of a potential well as a rate process with non- 
Markovian fluctuating rate, when the potential barrier 
exceeds ksT [54]. Such an escape dynamics is asymp- 
totically stretched exponential, but it tends gradually to 
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a single-exponential with the increase of the potential 
barrier. Then the rate fluctuations gradually vanish and 
for sufficiently high barriers (exceeding e.g. 12 kgT 
for a - 0.5, and about 9 ksT for a = 0.75) the es- 
cape dynamics is excellently approximated by the cele- 



brated non-Markovian rate expression [60. .86i |9U 19211 . 
as it was shown by us recently ifs^ll . However, even for 
such high potential barriers, when the escape dynamics 
out of a potential well becomes practically exponential, 
the diffusion in a periodic potential is not normal, but 
anomalously slow. This is because of viscoelastic sub- 
diffusion is based not on the anomaly of the residence 
time distributions (i.e. divergent, or extremally large 
mean residence time, like in the case of CTRW subdif- 
fusion and akin mechanisms), but on the long-time anti- 
correlations in the particle displacements (and velocity 
alternations). This fact in combination with ergodicity 
makes GLE subdiffusion physically much more appeal- 
ing scenario, especially in biological appUcations where 
it can be combined with the quenched disorder due to 
the medium's inhomogeneity (not necessarily leading 
alone to the emergence of subdiffusion). 

The subdiffusion in periodic potentials is described 
by some time-dependent a^ff which relaxes very slowly 
to a. The time pace of this relaxation depends very es- 
sential on the potential amplitude fsT], even if asymp- 
totically this subdiffusion is not sensitive to the pres- 
ence of potential |8lJ] because of the antipersistent, slug- 
gish nature of this subdiffusion finally wins over ther- 
mally assisted hops between potential wells. The latter 
ones can take place quite frequently and do not pro- 
vide a transport limiting step asymptotically (a huge 
difference with the CTRW approach!). This circum- 
stance makes asymptotic ratchet effect impossible in the 
limit of adiabatically slow driving with vanishing driv- 
ing frequency, O — > 0. A somewhat similar suppression 
of the response to adiabatically slow driving was ob- 
served also in the non-Markovian stochastic resonance 
[98^ ,99J . It seems to be a general feature of the non- 
Markovian dynamics with infinite memory which cul- 
minates in the death of linear (and not only !) response in 
the case of subdiffusive CTRW and the associated frac- 
tional Fokker-Planck dynamics |4ll-43|. However, for 
a finite Q the asymptotic anomalous current response 
of subdiffusive ratchet dynamics is present. In the low- 
est order of driving strength A it emerges as a nonlinear 
quadratic response to driving which violates the symme- 
try of thermal detailed balance in the absence of driving, 
for the unbiased potential. This is similar to the out of 
equilibrium physics of normal diffusion ratchets. The 
rectification current is optimized not only for a driving 
strength between two critical values Ai = 3nUo/i2L) 



and A2 -2A\ for the consider potential, but also versus 
Q. and t/o for the fixed temperature T . This is because 
of the rectification current vanishes in all these limits: 
n — > 0, n — > 00, ;/o — > 0, f/o — > The visual 
inspection of the single particle trajectories shows that 
the maximum of subvelocity versus Q corresponds to an 
optimal synchronization of the particle motions between 
potential wells and the potential tilts. This mechanism is 
somewhat similar to the current optimization for flash- 
ing, or pulsing normal diffusion ratchets |3] . However, 
the physics is different. In our case, the effect results 
because of the potential is not felt strongly by particles 
traveling large distances between the potential alterna- 
tions. Moreover, for sufficiently large frequencies parti- 
cles subdiffuse in the counter-intuitive direction, i.e. the 
anomalous current inversion occurs, cf. Figs. [3] El Like 
in the case of normal diffusion rocking ratchets |2], this 
effect does not have a simple intuitive explanation. 

To conclude, this work put forward a generalization 
of the pioneering contributions Jil [2l to the realm of 
subdiffusive Brownian ratchets in viscoelastic media 
rocked by a periodic force. The author is confident 
that the bulk of future research on subdiffusive Brow- 
nian ratchets is ahead because of their surprising and 
counter-intuitive features which call for experimental 
verification. 
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